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I. INTRODUCTION 


in ааз turbines the need to increase thermal efficiency 


f 


and specific power 15 a major concern. Опе way ШӘ τη 


Increase turbine inlet temperatures. AL pr 


№ 


sent, increasina 


sucn temperatures is limited by thermal fatigue strengths, 
Creep strengths, and melting points of the alloys used in 


turbine components. While the deveionment of improved alloys 


Nc. гаі Οἱ the salution, the development of efficient 


> = > =, + = -— . ` + Eo m ک‎ = 1 Е $ - 
соо! па svziems oi lade anj enga’ l] surfaces 16 just iz 
: UNE 1. vU E 4 т ων η : UE EE R πα 1 

important CRef. П ne, oH ας мвюс-іпп ап еттісісепіт === и 6 


ul 


configuration, the heat transfer distributians for the ga 
turbine components are needed. Because of the complex 


geometries and flows involved near blades and  endwall 


ul 


accurate convective heat transfer rate are meted σαι. E 


μὴ 
ІЛ 
о 


Flows through turbine passages are highly complex, 


containing numerous vortices, and secondary tlows. as well 


4 


as threae-dimensionai seperations. In recent studies, large 
portions of these complex flows have been identified and 
mapped out. The sketch in Figure 1, taken from L[Ref. nuc 
shows the flows which are thought to exist. As the inlet 
boundary layer approaches the blade, one portion forms cross 


ie meena the ather approaches the blade. Just in front of 


the blade a horseshoe vortex forms. nt the saddle point, it 
splits into a vortex on the suction side and a vortex on the 
pressure side. The pressure side vortex becomes the passage 


vortex, moving from the leading edge of one blade towards 


the low pressure side of the adjacent blade. The number of 
rotations of the passage vortex have been exaggerated in 
Figure 1 for clarity. Ordinarily t rotates Onl; ош ео 


twice through the passage. As the suction side vortex 


convects along the blade, it is eventually pushed away from 
the  endwal! ty the passage vortex from an adjacent blade. 
(Ref . 22 А smaller, corner vortex, rotates in an opposite 


direction to the passage vortex as verfied Бу Sieverding 


(ТБЗ 7 Get pa E 


B. Γι COG ING 

Film cooling is used to provide cool fluid between a 
surface and high temperature free stream gases to which it 
is exposed. Film cooling not only insulates the surface but 
acts as a heat sink for the hot free stream gases. The 
Overall effect of film cooling is to reduce the temperature 
of the developing boundary layer, which in turn reduces heat 
transfer to the surface. 


In our studies, knowing heat flux, q”, h is found using 


dEn I 957. 


The effect of film cooling will be indicated by the h 
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s PUL on for constant q'and Tig while varying Ty. When 
Ти = To + equation 1.2 indicates αἲξ 0, which may not Бе 
ας καν correct with film cooling. 


Another approach is one used by Goldstein, (Ref. 131: 


q =h Tw- Taw?) Γη 


Here 4” = 0 when Tw = Tawwith film cooling, by defintion- 


Without film cooling Trgm= Tay and the two equations are the 


same. Effectiveness of film cooling may then be defined 
Using 
n= Jaw -Tr 
Toc Тү CPAS) 


The [ἴπεθί. 1] method is not used for this study because of 
the difficulty in finding Taw for the test surface. 

The EHI s.s τ, ШІТРЕБЕНЗ protect a 
surface is dependent on numerous factors. The most important 
of these are: the blowing rate, size and number of injection 
holes, the location of the injection system in relation to 
the affected surface, and the curvature of the turbine blade 
or surface from which the injectant issues. The curvature of 
the blade in relation to the blowing rate is very inportant 
in the effectiveness of the cooling system. At low blowing 
rates, film cooling is more effective near convex surfaces 
than near flat surfaces or concave surfaces. If the blowing 
rate 15 Increased, this behavior is reversed such that 


concave surfaces exhibit the best performance. 
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REPRODDCEb"^AT GOVERNMENT" EXPENSE 


The approaches given ES equations 1.1 anda 1. еее 
compared for constant property "лен by showing linear 
dependence of h/h on 0, where Ө is the nondimensional 
temperature of the film coolant. Equating equations 1.2 and 


1.3 gives 


һ-һ,(1-70) Ca) 


Dividing Бу Бо then yields. (Ref. 4] 


ea C179) 
ho ho 


G OTHER SIUDIES 

Numerous studies have been conducted on the effects of 
film cooling on heat transfer, effects of secondary flows on 
heat transfer, and more recently, on the effects of film 
cooling and secondary flows on heat transfer in a turbulent 
boundary layer. 

Studing the effects of secondary flows on heat transfer, 
Gaugler and Russell (1984) compared visualized secondary 
flow patterns Over Six vanes with heat transfer 
distributions. The authors found that the horseshoe vortex 
Causes а local peak in the heat БАЗЕ ТЕ. rate near the vane 
leading edge. Large peaks in Stanton number меге found 
downstream of the vanes but the authors concluded that this 
effect was related more to vane wake rather than the induced 


horseshoe vortex. [Ref. 5] 
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Studing the effects of film cooling and secondary flows 
оп heat transfer, Goldstein and Chen (1985) ВЕН ап 
experimental study on the influence of the endwall on film 
cooling of gas turbine blades using a single row of 
Imyection holes. The authors concluded that the convex side 


ої the blade was highly effected by the flow originating 


from near endwall.Here film coolant was swept away by the 


passage vortex: whereas the concave side was not 
Sami t teant iy effected Бу flows originating near the 
Επηπαὶἰ |, Softee ο ας επ ihe affects of fila cooling 


ку 


Юе en ic Refa 851, ІБеі. 71. and LEef. 81. 
Studing the effects of secondary flows on heat transfer, 
Eibeck and Eaton (1985) conducted a study of a single vortex 


embedded in a turbulent boundry layer over a constant heat 


flux flat plate. The authors found significant increases and 


= 
id 


ni 


ES. im lotal Stanton numbers, due to thinning of the 
toundry layer оп the  downwash side pf LS егте Tand 
thickening оп the upwash side of the vortex.  Spanwise heat 
transter was larger as the circulation of the embedded 
vortices increased. [Ret. 9] Other work on the effects otf 
embedded vortices is given in [Ref. 101 and [ποί, 111 
Ongdren (1981) performed a study on the effects of film 
zooling on the heat transfer on the endwall of a turbine 


ascade, where the effects of secondary flows and cooling 


injection rates were observed. The author found that 


measurements showed a good match with predications in the 
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zase of no injection. Heat transfer measured on the center 
streamline of the endwall showed Stet i tative agreement with 
that of the flat plate. Secondary flows, which were believed 
to cause non-uniform distributions of film coolant, resulted 
in variations іп heat transfer across the width of the 


endwall. CRef. 2] 


Din PRESENT ΓΣ; 


The object of the present study is to increase the 
understanding of the effects of a longitudinal vortex pn 
heat transfer in a film сосіеа boundary layer. Th=#ëse ertTec-+s 
ar2 important regarding turbine blade encwall heat vransrer. 


The experimental model was constructed in a series of steps. 
The first was the design and construction of the in jest ren 
system, including qualification testing to verify uniformity 
and controllability. as well as the calculation of discharge 
coefficients for ihe various flow conditions. The next step 
was the design and construction of the heat transfer test 
surface followed by qualification tests to verify uniform 
neat flux and an energy balance to identify and quantify the 
sources of heat loss. After completion of the experimental 
apparatus four types of tests were conducted: heat transfer 


data with developing boundary layer only. with boundary 


ni 


layer and injected film cooling, with boundary layer and 
embedded vortex and finally with boundary layer, film 


cooling, and embedded longitodinal vortex. 


АЕ МЕТРЕ Επ μας. 


p ΕΠΟΝ SYSTEM 

The injection system was designed and developed to 
provide film coolant at temperatures above ambient. The 
Zolani is ejected from a single row of injection holes into 
the boundary layer developing along the bottom wall of the 
wind tunnel. The diameters of the injection holes were 
scaled relative to boundary layer thickness to be similar to 
з turbine blade, with a O/d ratio ranging Iron- Oa (Пп 


DD The tree str 


m 
Ш 


m air is at ambient temperature and thus 
the direction of heat transfer is opposite that of а gas. 


turbine. The temperature range for this study has been kept 


ο 


tJ 


mall ΚΟ. = Sts) to minimize the effects oT variable 


ІЛ 


properties. 

The injection parameters m and Ü were scaled to resemble 
parameters near gas turbine blades where m ranges from 09.5 
to 1.0 and Ó ranges from 1.2 to 1.8. Due to the reversed 
cc og οί heat flow for our πηι apparatus the 
Tre ‘Tro ‘Tw ratio is 1.04 - 1.07:0.74 - 0.95:1.0 as compared 
to cL ου στ αι. ὁ του actual as turbines. 

ІІ Ὀδεσςιβτίσῃ 

Bar for the injection system originates in ап 
Ingersoll-Rand air compressor, (two stage, 130 BSI, ο το, 


mode ] number ТІТП), where it is then sent to three large 


storage tanks. As Shown in Figure 2, the air then flows 
through an adjustable pressure regulator, a cut off: valve, 
reinforced flexible tubing (2.54 ст, і іп, inside diameter), 
moisture  seperator, flow regulator, a Fisher and Fortor 
rotometer (full scale 9.345E-3 m3/sec, 19.8 TSOR ΤΕ- 
number 104А3569А), a diffuser, and finally into the injecticn 
heat exchanger and plenum chamber. The rotometer monitors 
the volumetric flow rate tor filim cool no παπα о 
the diffuser is shown (п Figure ss 


А photograph of the chamber Ts эрозия свое лш 


The chamber is monstrocted of 1.27 cm СІ ПТР ЕКЕ Е 
with outside dimensvens gt 0.305 Оо ee ys 
18 In). As shown cm Emo ee the internal structure 


consists of three thin metal plates 0351000200 Tete 1 
Г) — 3.08 cm (2 In) trom the notrom and рес πα 
ор at 3.08 cm intervals. Two silcon Yebber heaters. у ш 

О. ея м τα е MEE 120 volt. аге seperately placed over 
the bottom two plates. The heaters are controlled throuah a 
rowerstat variable autotransformer (Lype 1955). ШЕ rap 
surface contains 13 plexiglass injection tubes each 8 cm 
long with an inside diameter of 0.99 cm (3/8 in) with a ΙΙ 
ratio of 8.42. The 13 injection holes are inclined at an 
anale of 30 deg., with 3 diameter spanwise spacing between 


center пее where the middle tube is set on the centerline 


Ot the гезиги аар 


HJ 
۲.3 


Air enters the chamber through the diffuser section 
and is then directed over the two heating surfaces where the 
air can be heated from ambient temperature up to BO °C by 
controllina the input voltage to the heaters through a 
Variable autotransformer. 

Three pressure taps, positioned at the center of the 
front and two side faces, are used to measure Рос = Pg. 
Three 0.254 mm (0.010 in) diameter copper-constantan wire 


tnermocouples with welded junctions are place at different 


lazations inside ihe chamber to measure the uniformity 


1-1 


Тор in the plenum. 
ERN ας Гг арго and Performance 

The uniformity of the plenum chamber pressure, Ре, 
was found to be satisfactory over the range of injection 
conditions with typical differences of approximatly 1% іп 
the spanwise direction and a maximum of "οτί. TOF 
only one case at a lov flow rate of 0.327E-4 m9 /sec. 

The plenum produces a reservoir of air at ап 
elevated temperature and pressure, which is near stagnation 
conditions. The temperature at the nozzle exit, Tre; is not 
equal to the temperature in the plenum chamber, Тор, due to 
conduction heat loss through the nozzle surfaces to the 
surrounding ambient air. l. is necessary ο know the 
relation between Тс and Top because injection parameters 
are estimated at exits of the jets and plenum chamber 


temperature is more convenient to measure. A test was 


or 
aL 


conducted to determine the relationship between these two 
temperatures over the flow rates ranging from 0.327Е-4 to 
O.701E-2 m/sec. Results are shown in Figure 6. The data 
trend shown in the figure was the same regardless of flow 
rate, and an average of exit and plenum temperatures меге 
used to determine each data point. Exit temperature shows a 
near linear dependence on plenum chamber temperature. A 
curve of the form Toc = ο ος where G = 1.455 апа В = 
0.858, was fitted to this data so that Tre could Бе 
estimated from plenum chamber temperatures in subsequent 
tests. 

The recovery temperature at the nozzle ех!ї, Тү, 15 
given by 

2 
lC 2571) 

where QŒ is the recovery factor, having typical values of 0.4 
to 0.9. QŒ represents the percent of dynamic temperature 


which is not lost to viscous dissipation, where Q is def ined 


by 


(252) 





The total nozzle exit temperature, Toc, may be expressed 
USING 


ué 


Toc = Te + `29сСр 
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Because of the low velocities employed, "(and neglible 
viscous disipation), Toc = Тіс within a fraction of a 
degree. 

Measuring parameters Ро: о. УЕ and Тос апа 
knowing the area, А, normal to the flow, of the injection 
holes, designed to be 9.2633E-4 m2, the coolant velocity Ue 


is given by 


За 
Uc? A (2.4) 
The static density ıs estimated using 
دم‎ в) 
ς R Тс 


Шввв (10%, mc. is now the product of Uc and pc. To calculate 


the isentropic mass flow, Юс; апа Uci are found using 


. .Poo 
Pci ^ RTop ο ο) 
and 
шен 2 Рос Poo) 2. 
Uci = Pci PL 


These two equations are derived from the equation for 


compressible flows: 


27 1 


2, 2 
Poc 4 7 Poc : ; 


Discharge coefficients, Cg, are then estimated using 





d (OcUc)j а. 


hF.) 
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The variation of discharge coefficients with 
volumetric flow rate is shown in Figure 7. The dependence is 
near linear, where Cg decreases slightly with Tc at a айық 
flow rate. Discharge coefficients range from 0.581 to 0.730, 
wh ich is consistent with other workers’ results (πεί. ο]. 
When plotted verses the Reynolds number as shown in Figure 
В, the discharge coefficients at ditferent temperatures 
collapse more closely together than in Εὐπῦτε πα. Susp» 
behavior indicates satisfactory injection system performance 
[Ref. 4]. Numerical tabulations of the data are found in 
Τα τες 152 ὃν αν D а 

Curing qualifizations tests small hysteresis was 
indicated by the plenum chamber pressure as the fiov 
inzreased and decreased. With increasing flow rate plenum 
chamber pressures relative to ambient pressures, were 2 to 
10% lawer than when flow rates decreased. These data were 
taken over time intervals long enough to ensure steady state 
conditions were obtained. This result was probably due to 
fluid Circuit effects Gi БЕ ТГ ШЕРІ ΝΠ, 
chamber, specifically. behavior similar to that клише ET 
inductances and capacitances in electrico circuits. This 
hysteris was probably the result of a velocity dependent 
Separation occuring in the diffuser of the injection system. 


However, in 


μὴ 


pite of the hysteresis, data are self 


consistent (even though obtained as the flow rate both 


decreased and increased) with minimal scatter: Са varations 
trom hysteresis are no greater than 1 or 2%. 


The heat transfer surface was designed and developed toa 
mous —soenstamnmt heat flux over its area. The average 
Surface temperature may be adjusted and maintained from 


ambient up to 60 С. The device was constructed so that the 


Upward facing part is exposed ta the wind tunnel airstream, 


with minimal near unes bwecondwsliilon тст пеш ПЕЕ а 502 
SGeneath t^s test μμ σα. dj crudos c ES has bazan 
instrumented to measure Surface temperatures with 


thermocouples place just beneath the surface, and a film otf 
ШІНШІ сгусізіс sprayed on the top surface. 
ШІ есггірііоп 

ñ photograph of the heat iranster surtace is shown 
in Figure 9. The design is based on a Similar surface used 
by the University of Minnesota L[Ref. I131 ang σεις ы Ла ГТ 
consists of a thin stainless steel foil τ. 3O2 TOI 
Ш. O. 127 mm x 1.194 m x 0.467 ms painted flat black 
ШЕТІ ος ἰπαιοῖν S layers of liquid crystals. Attached to 
the underside of the foil are 88 copper-constantan, 0.254 mm 
diameter, thermocouples integrated into grooves cut into a 


EE hop οἱ 0.25 mm (10 mil) thick double-sided tape 


(manufactored by the 3M Company). The grooves are then 
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filled with RTV. А silicon rubber heater, 1.0 mm x 1.143m x 
0.457 m, 120V/100U, is attached to the тара wit a ENTER  ;; 
epoxy which was applied to both the tape and heater 
surfaces. The heater is then backed by а 12.7 mm (1/2 DITE 


thick lexan sheet, followed by 235.4 mm of foam insulation, 


ΓΕ 


2.85 mn thick eStyroftoum ana ӘНЕ ΠΕ Ἐν, ή tu 
balsa Wood, as shown in Figure 10. 

Around the edges of the foil grease was inserted 
between tha foil and plexiqlass frame to allow for thermal 


re also 


fu 


expansion oi the top surface. Twos nin ШЕКЕ 
attached to tne back corners of the roi Pr s; jp o, 
through the floor “or Stile ward etum sl E lo We Veins 
attached to each wire to add tension to the foil and help 
maintain a flat surface during thermal expansion. Additional 
vertical movement of the foil surface above the bottom wall 
of the wind tunnel occurs due to thermal expansion during 
heating. The surface is maintained level by adjusting screws 
in the plexiglass frame supporting the heat transfer surface 
from below. During heat transfer tests, the top surtace of 
the foil remained remarkably flat and smooth with minimal 
surtace irregularities. 

Thermocouples are placed on the surftace as shown in 
Figure il. In each row, thermocouples are 2.54 cm apart. The 
surface temperature is controlled by adjusting input voltage 
to the heater using Е Standard Electrical Product Co. 


variac, type 3O000B. Uith this type of heat transfer surface, 


"ο 
ra 


audition of more thermocouples for improved spatial 


ο ποπ of Surtace temperature is difficult because of 


cold spots and paths for leakage caused by large numbers of 


thermocouple wires between the foil and heater. 
<. Qualification and Ferformance 


- 


Preliminary qualification tests were made using ап 


ο version οἱ thse heat transfer surface in order to 
understand its behavior and performance characteristics. 
Improvements in design and construction led to the heat 
transter surface used to obtain the final results. 
ЕЕ ИТ тат гоп results from the final surface, using 


temperature sensitive liquid crystals, showed more uniform 
temperatures. giving evidence of more uniform heat flux, 
than results from the early version. 

To test the early version of the heat transfer 
surface, a Hughes FProbeye Thermal Video System series 4600, 
РЕСП ot a infrared and video camera with display 
screen, was used to measure the surtace temperatures with 
heat Tramster surface outside the wind tunnel. With this 
system temperature variations as small as 1 °C can be 
measured. The surface was observed Raden four operating 
renditions: natural convection with a surface temperature of 


m 


approximately 33 os forced convection with ап initial 


қ Р | E o 
surface temperature of approximately 33 С, пата! 


: e 
convection with a surface temperature of approximately 45 C. 


and On ed convection with ап initial temperature ої 
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approximately 40 °С. Results of these tests are shown in 
Figures оис sita 15 respectivily. The test were 
undertaken to measure temperature differences in order to 
provide some verification of uniform heat flux under 
different operating conditions. 
a. Preliminary Tests 

The natural convection conditions shown in 
Figuras 12 and 14 indicate that the surface temperature 
distribution (susmi ppo li Im = 3 С. ñ low temperature 
region is located near the leading edge on the left side, 
caused by the larger number of thermocouple wires located 
near the forward section of the test surface. This problem 
was eliminated in the final heat transfer surface design. 

Figures 13 and 15 show the surface isobars for 
the forced convection tests to be spanwise uniform midway 
between the two sides along most of the length of the test 
surface. Such behavior evidences a locally two-dimensional 
turbulent boundary layer with a constant heat flux at tne 
surface. Moving air was provided by a Таһов Preduct= sine. 
fan (model number PF-O0S-1) with a 22.25 =neblasewe di anette. 
which produced а free sity eam velocity estimated to be 
between 72 and 4 m/s. Away fram the flow produced by the fan, 
near the sides of the test surface, temperature variations 
are highly non-uniform as would Бе expected. Comparing 
results in Figures 12 and 14 with those in Figures 13 ала 


15, the small cold Spot present with паәшшга)! EOE IE ш: 


= 0 


the leading edge almost disapears with forced convection, 


Probably а5 a result of the higher heat iransfer 
si I ents involved, 

А ене па cool spot is located near the 
το πο edge and is present with all surface conditions as 


τη Figures 12, 14, and 15. The variation is caused by 
a small non heated portion of the silicon heater containing 
aF ir mosv.at. n the final design, this heater was replaced 
with one without a thermostat to avoid formation of a local 
ENG Spot. 
En Energy ba lange 

An energy talance was performed to determine 
heat loss by conduction оп the heat transfer surface used 
to obtain final results. During the energy balance, heat 


loss. Бу radiation and convection were prevented since the 


metal foil surface, ordinarily exposed to convection in the 
ο. tumnel9 was covered three layers of, 25,4 пт thick, 
foam Шеш Om. τον the “energy balance and for all wind 


tunnel testing, foam insulation was also placed around the 
sides of the test surface located below the wind tunnel 
convection surface. To estimate heat loss through the 
СЕЕ Ва оп "placed on top of the foi! surface, the one- 
dimensional, linear farm of Fourier’s conduction equation 


wae employed 


а = КА =— (2108 
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REPRODUCED AT GOVERNMENT EXPENSE 
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where k is .04 W/m K, A is 0.4897 m2, X is 0.0254 m, and AT 
is the temperature drop in the X direction. Heat conduction 
through the bottom and sides of the heat transfer device are 


then given by 


acond ° VI - qç ІІ 
where VI is the power into the test plate, and q is the 
conduction loss through the top insulation. Tests to 


determine Qeong were made at five power levels 12.00, 14.20, 
15.52, 20.44. апа 20.45 Watts. These levels were choosen to 
te of the same magnitude as conduction losses experienced 
under normal test conditions with convective neat transter 
from the top of the plate. 

acond is plotted verses T, - T amb in Figure 


17. From these results an emperical equation was generated 


Чсопа 9.93 * 1.45 АТ - 0.051AT? + 0.00068AT3 (2. 125 


Where, here AT is Ty = T amb. This equation is only να 
over a range of Tw, - T amb from 8 to 30^C. Шһеп exposed to 
convection in the wind tunnel, conduction losses are only 
1.5 to 2.5% of the total power, therefore even a 25% error 
in the estimate of conduction losses will cause less than a 
1/2% change in the estimate of heat transter by convection. 
In order to have additional verification of the 
conduction model given by equation (2.12), conduction losses 


were estimated beneath the surface, through the plexiglass 


sides, as well as from the foil to the plexiglass frame. 
Numerous thermocouples were placed at key positions between 
the layers of material and the test plate to make these 
estimates, as shown in Figure 16. Temperatures from these 
thermocouples were then employed in the one-dimensional form 
of  Fourier's conduction equation. The estimates were made 
for the same power levels mentioned earlier. Results аге 
listed in Table 6. The total of these estimates showed rough 
agreement with equation (2.12). At low power inputs, less 


than 16 Waits, the mode) slightly under estimated losses, 


while at higher power levels, greater than 18 Watts, the 
model slightly over estimated conduction losses. The 
differences are probably due to the multi-dimensional 


effects at the corners and sides of the foil surface, which 
could not be accounted for using a one-dimesional conduction 
equation. 

Radiation losses were estimated from the 


following relationships (Ref. 15] 


Аа 
аа Сі. еі j VIE T 
А, А ij €jAj 


and 
N 
Grad =2_ij (27140 
1=1 


The view factors, Fij ; for the top wall and each side wall 


33 


REPRODUCED AT GOVERNMENT EXPENSE 


KVIISIssOoI€KCSCIK — EU, to EA C EN. cc ee ee 
- = оз ое = = 


— 


were estimated to be 0.66 апа 0.17 respectivily, where 
N 
2 Fl 1213) 
j=1 


From these equations, radiation losses for an average plate 
temperature of 40°C was found to be 53 Watts, approximately 
8.5% of the total power into the transfer surface. 
с. Contact Resistance 

In the present study, thermocouples are attached 
to the stainless steel foil using the 0.25 mm double-sided 
tape and RTV silicon rubber epoxy. Welding was not used due 
to the thinness of the foil. Consequently, a thermal contact 
resistance, 1/(НсА), is present between the back of the foil 
surface and the thermocouples. This contact resistance along 
with the thermal conductivity, К, of the meta! foil will 
Cause the thermocouples to measure a higher temperature than 
actually found on the test surface, where this difference, 


AT, is given by 


= 1 АХ 
atsa (ir М 4%) 215) 
Contact resistance is highly dependent on the contact 


pressure as well as the area of the surfaces in contact. Due 
to the method of construction of the test surface, contact 
pressure can only be estimated while the area of the 
thermocouples in contact with the surface can vary in the 
shape and size of each bead. These unknown properties 


precluded an emperical estimate of the value of the contact 
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resistance. Two different methods were used to determine the 
average contact resistance experimentally: (1) with 
conductive heat transfer conditions and thermocouples above 
and below the foil, noting that the contact resistance is 
twice аа large due to the presence of two layers of 


thermocouples, and (2) by direct measurement of the drop in 


temperature with convective heat transfer. The. second method 


involved comparing the surface temperatures indicated by 
liquid crystals to ones measured with thermocouples beneath 
the surracese. From these two methods the contact resistance 
plus the conductivity; resistance in the plate was estimated 


he ος I al The Eme UE уаз» used for al} 


i 
[1 


thermocouples. Because this is an average value over the 


entire plate, the actual resistance for individual 
thermocouples may vary from this value, causing small 
deviations in the spanwise heat transfer coefficients and 


STantonmenunbers . As μου] be expected these variations wera 
unaffected Εν ГОУ condit i Ons. ws wil De seen, small 
variations due to contact resistance are eliminated Бу 
presenting results as Stanton number ratios for local 
вов tions. 
2 МЕ arl ver Deslan 

The final heat transfer surface design is the third 
iteration. The first one had 146 thermocouples held in place 
by two layers of tape. Due to the large number of 


thermocouples and associated wires, the surface did not 
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exhibit uniform heat flux. The thermal video system 
indicated that isobars were nonuniform in the spanwise 
direction with many cool spots above the locations of large 
numbers of thermocouple wires. The cool spots were observed 
under both natural and forced convection.: The оспа UN 


A} 


ІЛ 


o had some small dents. incurred “resi рр Иб CO Бекжан 
grade of stainless steel foil, packaged as a single sheet 
rather than a coil, was used in later designs to minimize 
surface deformities. 

Of the three designe. the final one was believed to 
nroduce the most uniform surface heat flux because (1) the 
surface was the smoothest, (2) fewer thermocouples were 
used. attached with tndusrial strength double backed liner 
instead of 2 layers at tape, (3) thermocouple wires were 
spread uniformly behind the foil instead of being кыр 2 in 
groups, and (4) slots were caretully cut into the tape for 
placement of thermocouple wires to prevent and raised 
surfaces between the foil and heater. A plexiglass frame was 
also added to provide better support from the sides in 
addition to a lexan sheet which was added to supply a more 
rigid support from below. The lexan helped prevent sagging 
and to maintain continous contact between the foil, tape, 


and heater. 
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C. TEMPERATURE MEASUREMENT 

ATI thermocouples used in this study were type-T 
opper- zonstantan), 0.254 mm diameter, manufactured Бу 
Omega Engineering. Two thermocouples (freestream probe and 
temperature profile probe) were calibrated USING a 
temperature regulated bath consisting of liquid nitrogen and 
electric ws. and a known platinum resistance reference 
0501 Cis The calibration was performed over the 
temperature range of 18 - 60 °C. The millivolti reading of 
the thermocouples verses, the reference temperature was used 
to Qenerate a fourth degree polynomial equation using a 


least squares curve fitting technique. 


T-26.573E - 1.937 E? * 0.998E? - 0.26 1E ЕП 


where 


Е = millivolts x 1000 СЕТЕ 


The same polynomial was used for the thermocouples used оп 
the heat transfer test plate as al] ео les indicated 
very similar behavior at different temperatures. 

Temperature profiles were performed using one of the 
calibrated thermocouples mounted on a traversing mechanism. 
The traversing mechanism was driven by a simple threaded 
shaft (18 threads per inch). This threaded shaft was used to 


move a traveling block attached to the mounting stem of the 


oa 


thermacouple probe, such that shaft оао p Sh ГЕ 
converted to vertical probe position with respect to the 


wall. 


D. WIND TUNNEL 
The wind tunnel picivred in Figure 168, Guilt Gy er 
was .desiganed to provide a flow field from the nozzle with 
uniform velocity and low turbulence intensity. P is 
designated the NPS Shear Layer Research Facility (SLRF). The 
tunnel has numerous pressure taps along both side walls and 
Tontains many instrument ports along the top wall to measure 
tT TOW σπα Εμμ ας. 
l. Description 
The SLRF, shown in Figures 20. is а wind 
tunnel of the open circuit blower type with fan upstream апа 
air entering the fan inlet from the surrounding E The 
air speed through the test section can be adjusted from 5 to 
40 m/s. The tunnel frame has leveling screws to adjust the 
centerline of the tunnel Та а “есіліп ган position. Тһе 
discharge part of the tan slips into the inlet end of the 
wide-angle diffuser with a 1.56 mm (1/6 in) clearance all the 
way around to isolate vibrations from the fan to the wind 
tunnel body. The diffuser SELTVen contains a filter pack and 
a nozzle leading to the test section. The test section is a 
rectangular duet, 37048 m (10 ft) Vong, апа ο ου ο πι Uc MEM 


wide, as shown in Figure 20. The top wall 18 а continous 
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panel fabricated from 4.76 mm thick Lexan sheet with 
neoprene channel seals at the edges. The ceiling panel 
height and shape may be changed to permit adjustment of 
Е pressure alang the length af the test section. A test 
section height variation of 0.1524 to 0.3048 meters (5 to 12 
ВЕ in be oetained. The floor of the test section 


Éronsisests of. 0.6096 meters (2 ft) and 1.2192 meters (4 TLD 


D 


sections which are removable and replaceable. Each floor 


section has "O" ring seals at their seams. Removable side 
ЕЕ allow Όσο accessibility to the {ast section, 
Profils measurements may be made using probes inserted 
through the top and bottom walls.tRef. 15] 

à schematic of the test section components used in 


the present study are shown in Figure 21. Ап unheated 
starting length of 1.10 m exists upsteam of the heated test 
surface. The injection nozzles are located 1.08 m downstream 


t 


2 


ої the boundry layer trip and 0.02 m upstream of the te 


n 


surface. The leading edges of the vortex generators are 


"ο 


placed 0.47 m downstream of the boundry layer trip. А 
schematic of vortex generator geometry is Shown ОП halen е 
r. 
2. Qualification and Performance 
Extensive Qualification test of the Shear Layer 
Research Facility were conducted By  Ligran! (Ref. [d 


Results show that the variation of total pressure at the 


exit plane of the nozzle is less than 0.4% at 2& m/s and 34 


m/s. Mean velocity varies less than 0.7% for ШЕ same mean 
freestream Speeds. From five-hole pressure probe 
measurements, the velocity angle deviation is nowhere 
greater than about 0.6 degrees at the nozzle exit plane. 

Profile measurements of the mean velocity and 
longitudinal turbulence intensity in the turbulent boundary 
layer developing at 20 m/s indicate normal, spanwise uniform 
ο. For this qualification test, and all results which 
τε] ν᾽ the boundary layer was tripped near the exit af the 
nozzle with a 1.9 mm high strip of tape ας πι 
measurements along the test section surface ai ihe nozzle 
exit were uniform within 0.55 indicating spanwise unirorm 
akin Trie Troms 

Freestream turbulence intensity was measured to te 
0.00083 (8.3 one - hundreths of опе percent) at ΟΤΙ 


increasing to 0.00095 at 30 m/s. 


Е. Шата ОсОО CTER 

The data acquisition system was designed to rapidily 
measure thermocouple yoltages and convert them to 
temperatures in deg C. Using these temperatures, along with 
user supplied information on ambient conditions, freesiream 
conditions, and power input. the system calculates free 
stream velocity, density, local heat transfer coefficients 
and Stanton numbers, and spanwise averages of heat transfér 


Goefiitients and Stanton nuntpen=_ 
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1|. Hardware 
The data acquisition system is comprised of an HP- 
SJB computer with its associated memory module and interface 
BEES. Storage is provided by an internal cassette 
storage drive. Voltages from the type-T (copper-constantan) 


thermocouple 


Л 


ате read by an He se to ih Data 
Acquisition/Contral Unit with a HP-3498A Extender. The unit 
communicates with the computer through a ΠΕ ο τῇ 
Interface. Software, described below, converts the voltage 
inputs into temperatures. The HP-S5B is used to "calculate, 


store. display. and print desired information. 
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 ireemiprosSrams, ‘ЭТОТ, 500815. and TPROF . were 
developed to be used under various flow conditions. The 
three programs are listed in Appendix D. 
SCENE DAT 
SIDS, developed Ly Ligrani, ӨКГІс; and. Joseph: 
prompts the user for current (amps) and voltage ta the 
heater to calculate the power (Watts) into the heat transter 
surface. The program continues by prompting the user ror 
stagnation pressure (in HO) of freestream and ambient 
pressure (іп На). The computer than reads the thermocoupie 
voltages, converting them into temperatures and storing them 
in a matrix. After all temperatures have been calculated the 


freestream density and velocity are calculated in SI units. 


ul 


Er omes point, the program accounts for conduction and 
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radiation losses. 
calculating local 


coeffwicients. and 


given based on the 


these calculations 
attached printer an 
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Stanton numbers. 


and contact resistance. . Itir oni ince REN 


and spanwise averaged heat 


downstream distance. 


transfer 


Reynolds number are alsa 


Upon completetion of 


the program prints the data оп ап 


d stores them on a disk file. 


is A modi ication: 


This program runs in 


for modifications to use 


baseline measurements Го оа ево 
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program also 
zoolina 


сЕ Гс епски ЕЕ 


injection 


contains an option 


parameters 


iy ratios. massflux 


ratios. and blowing rate. 


c.c TIL 


T PEDIS 


temperatures from thermocouple probe output, 


was traversed in the turbulent boundary 


including 


Бу Joseph, to the 
the same manner as 
previously stored 
type relationships. 
to "calculate 2777 60 


discharge 


ratios, momentum tlux 


was written Бу  Ligrani to Calculate 


as the probe 


layer, above the 


heated test surface. As the program begins, local flow, wall 


heat TUS and wal 
then calculates a 
velot itni The 


| temperatures are read in. The program 


friction 


program then enters 


temperature and a 


f ieia 


a loop where | 


acquisitions a voltage from the thermocouple, and then 


σα) να Ооп тесш е 


temperatures from the volatage usina probe 


. This loop is repeated as the probe is 
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moved to new locations, where probe location is read in and 
accounted for. Ad the last profile point (the probe is 
moved from the freestream to the wall), the program prints 
out relevant parameters including thermal boundary layer 


thicknes 


as well as dimensional and normalized profile 


d 
(12 
rr 
ni 
+ 
г 
3 
i} 
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III. EXPERIMENTAL КЕС Е 





The study was conducted in four partis. Те тте. 
consisted of establishing baseline Stanton numbers and heat 
transfer coefficients for a turbulent boundary layer over а 
constant heat flox test  suriace at Ug =10 m/s with an 
unheated starting length of 1.10 meters. The second part 
consisted Of documenting the effects of various size 
vortices over the same surface using the same freestream 


4 


ον, πα πω. тпе зе втер Involved measurements in 


Hi 


surbutlant Boundary layar with Trim тосо pon Rir e Ы 


4 


rates. The final series of tests were conducted to study the 
effects of an embedded vortex on heat transfer in the t dm 


zooled türbGulent πο υπ σος С” 


А. ΕΘΕΤΕ ΙΕ ΟΕ Εν» 
Heat transfer measurements were performed and recorded 


τος three different freestream velocities over the test 


A IE MT "ο. 19 and 20 m/s. Numerical гава с тБа OJ heat 
transfer coefficients, Stanton numbers at different 
thermocouple positions in addition to spanwise averaged 


- 


Stanton numbers By thermocouple row are listed in Appendix 
ES. 
The spanwise heat transfer coefficients of the test 


surface at 10 and 20 m/s are shown in Figures 23 and 24 


respectivily. Есери геги еси ihe spanwise uniformity is 
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very good with maximum variations of 10% (from the average 
for a given row) at 10 m/s. These small variations are due 
to differences in contact resistance between different 
thermocouples. As shown in the two figures the qualitative 
form of these variations is the same at different flow 
conditions. Larger spanwise variations for row 1 appear to 
be due to multi-dimensional conduction losses near the 
leading edqe of the test surface in addition to contact 
resistance. 

Spanwise averaged Stanton numbers for 10, 15, and 20 m/s 


are plotted as functions of Reynolds number in Figures 25, 


DS. and 27. These data show agreement with the empirical 
equation for turbulent boundary layers at constant free 
stream velocity along a flat plate, with constant heat flux 


and an unheated starting length of 1.10 m. CRef. 183 


Τμ. 
St,Pr = 0.030Re, Ë (=) | (3.1) 

The maximum variation between the emperical equation апа 
measured data is approximately 3%, at 15 m/s, with excellent 
agreement at 10 and 20 m/s. At a given Reynolds number and 
test plate location, Stanton numbers may be repeated within 
a few percent provided thermal equilibrium of the wind 
tunnel and heated test surtace have been achieved. 

Mean temperature profiles were measured at Х=1.44 and 
Icy meters. These аге shown in Figures 28 “and 29, 


indicating expected behavior. At 1.85 meters, profiles show 


REPRODUCED AT GOVERNMENT EXPENSE 


excellent spanwise uniformity at three different span 
locations, which indicates that the mean flow field is two- 
dimensional 

Figure 30 shows mean temperature data plotted in non- 
dimensional wall coordinates Y+ and T+. When compared to the 


empirical law of the wal) CrRef. 183 


T*=2.195InY* +13.2Pr - 5.66 (3.2) 


the data show agreement for 100 < Y+ <200. At larger Y+, 
data 28 behavior typical of outer Срег νους ὃν, thermal 
Ms developing with an unheated starting length. Mean 
temperature profile data deviate from equation (3.2) for Y+ 
£100 due to probe spatial resolution effects and significant 


heat transfer from the heated wall to the probe. 


В. SINGLE VORTEX 
Four different size vortex generators, Figure 22, were 


individually positioned 0.479 m downstream from the boundary 


layer tripo ΘΠ results were obtained at a free stream 
velocity of 10 m/s. The heat transfer are given in terms of 
St/Stgfor thermocouple position. The measured data for each 


test is found in Appendix G. 
The St/Stg ratios, as a function of 2 position, where 


2=0 is the centerline, are shown for each vortex in Figures 
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Vortex generator 81, smallest, was positioned with its 
mx EM Lett ot the wind tunnel centerline. The effects 


Meme vortex. are evident in the Stanton number ratios. The 


ratio increases near the downwash side of the vortex to 


ш 


maximum οὐ 1.19 and decreases over the upwash side to a 
minimum ої 0.975. The downwash side of the vortex appears to 
ce thinning the boundary layer, thus increasing the 
localized heat transfer, whiie the upwash side appears to be 


thickening the Boundary layer, decreasing the localized heat 


transfer. The baseline Stanton number, Si, data used to non- 
fegemearom@a!ize results in Figures 28-35; and all subsequent 
ΕΕ Wee Obta:inedg áit conditions which were not ait ёхагпштї 


thermal equilibrium. A repeat of these data give St, values 


. 
. 


1 


f 


а Там percent higher, which gives slightly lower 5%/5% 


SU Sits values at 


+ 
fp 
т 
5 
ul 
Ci 


Бететріпа to Figures 28-3 
eee ons away from the vortex are thus closer to 1.0 than 


Epnwisated in the figure 


ІЛ 


Vortex generator #2 was positioned 41.797 cm left ot the 
wind tunnel Senter bine. The е үс сы пел їгапе=?ег are 


πι σα the == 


ш 


me as with #1 only with ап increase in the 
ии от the St/St ratio to a high of 1.204 and a low 
of 0.949. The overall effect in the spanwise averaged 
Stanton number is nearly ihe same as with vortex #1: both 


are approximatley 5% greater than the baseline data without 


ar tex, 
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των ον. 43.9. 9. — 


COOL — £.* 45.9. лғ.» алда, лла . 
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Vortex generator #3 was positione 8.08 cm left of the 
centerline. Again, the results are very similar to those 
from generator Hi,with a high of 1.274 and a low of 0.951 in 


Stanton number ratios. Тһе spanwise average of Stanton 


numbers are nearly identical to those of vortex #1 and R2 
The major difference noted with vortex #3 is the smearing ot 
Stanton number ratios on the USWaeh siege ot спе vorte see 


positions 2.24 and 5.08 cm. This effect may suggest the 


presence of a small secondary vortex, 


κο 
wl 
a 
Uf 
P 
=) 
m 
—. 
«4 
Г) 
+ 


Vortex 


ΠΝ 


eneratar #4. SOS) vl ores 


centerline, soreduces a vortex so large that it daminates the 


entire measured fiow field of the surtace. The diameter of 
the vortex | appears to range from 10 to 15 “m. while the 
spanwise averaged Stanton numbers show a marked increase 


from the baselina of approximately 10 -14%. The eftects of a 


ul 


secondary vortex have become more visible and can Бе seen 


li 


0 (262 Ein. Due to the | 


—h 
о 
3 
[ 

H 
e 
co οἱ 


roe size ang effects ої 


vortex #4 it w not used in any subsequent tests. 


{11 
iri 


These results, especially those found in Figures 30 and 
3i compare very favorably with those of Ertech and Eaton, 
Figures 26 and 37. Figure 35 shows the spanwise Stanton 
number ratio wiih a vortex generator Z2 cm high and a ZÓ deg 
angle of attack. Figure 37 shows the spanwise Stanton number 
ratio with a vortex generator 3 cm high and a 12 deg angle 


о ак аск οσο MEN 
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The strength of the individual vortices can be 


quantified by their respective circulations, where 


circulation = 2.5HU 0-10 =) 


mores Study the circulation of vortices #1, #2, ЖЗ; and 
#4 were determined to be approximately 0.045, 0.075, 0.125, 


πο ο 175 m /s respectivily. (Ref. 191 


ШІ ЕТІП COOLING 

Heat transfer measurements were períormed and recorded 
Pens with film cooling at O= 1.543, 1.485, and 1.474 
Е ото ratios, 0.68, 0.55, and 0.54 respectivily. 
Numerical tabulations of local heat transfer coefficients, 
local St/Stg ratios, and spanwise averaged Stanton numbers 
are listed in Appendix F. Тһе Stanton number results for 
each of the three blowing ratios were stored in a data file 
to be used for comparison with those measured under the 
па впПсвЕ of toth a vortex and film cooling. 

When comparing this film cooling data with results from 
injection system qualification, measured discharge 
coefficients, were afew percent lower for the same flow 
rates. This was due to two seperate effects: (1) no external 
flow existed with qualification tests whereas in later 
tests, jets were subject to external flow as the exited, and 
(2) when the injection system is connected to the wind 


tunnel test section floor, small steps were present at seam 
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ПОСЕ σος in injection tubes. The slightly lower discharge 
coefficients, measured during wind tunnel testing, had no 
significant effect оп flow behavior in the film cooled 
boundary layer or on the accuracy of measured injection flow 
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show trends 
turbulent 

пјес es 
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blowing Fat bo set 

without 


η ια Ες 


apparent near wall 


deficits which extend to greater Y 


St/Ste ratios are -plotted E 
in Figure 32. Here a Reynolds 


to απ Αη ratio s s 


zcoOrrespomnj3s ооа raino som 


© edae of 


from the down Lream 


ion hole diameter d=0.95 cm. AS 


ihe overall heat transfer rate 


when compared with that of 


due 


the leading edge af 


consistent with other 


boundary layers cooled using а 


үе. c 


measured at 10 m/s with film 


ο ασ ang б -1.543. These are 


film cooling for the same 


The profiles with frim 


erences fram the freestream 


qradients and 


than withaut 


соо! тля. Such behavior evidences some effects produced 


by film cooling: thicker thermal boundary layers and lower 
heat transfer from the wall to the freestream. 
РЕМ ЕНЕГЕ VORTEX AND FILM COOLING 


Ihe final pnase of research consisted of measuring tha 


=+; 
iti 
do 


ps 


—, 


it 


Z e a vor isx neat transtiter in 


single 


turbulent boundary layer. Measurements were performed using 


4- 


E - enseratonr ne plated at three  dirferent locations, 
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Ea 


Were 


mi 


tu 
P^ 


Il 


22 


iÙ 
Lt 


in a 


rt 


w 


ΓΙ 
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h and m values. 
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Ап 


correspond to a boundary 


circular symbols 


undisturbed turbulent 


ле as Hume tron от 


а 


-4.79 cm, arae shawn п 


these figures, SUS 


boundary 
with 


layer 


shaw results for a Ec HAY 


layer with film cooling and an embedded vortex. Ав shown іп 
the figures, except for local hot spots, overall heat 
transfer rates are lower than in an undisturbed turbulent 


boundary. Localized heat transfer rates. greater than 
St/St=1.0 occur near the vortex downwash side. The increase 
begins at approximately A=1.9 meters and continues down the 
length of the test surface. At A=1.457 meters the maximum 


the 


Ut 


where S&St/St. =I τσ. ne teen interesting feature i 
overall drop in the Stanton number ratio, гапа па ГОС 
11%, on the upwash side of the vorter.: This етет вере 
not anly in the downstream direction along the length of the 
surface, but also in the spanwise direction. 

Figure 30 shows the effect of an embedded vortex on 
Stanton numbers in a film cooled tore Eoeunuum ye E 


figure it is apparent that there 
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15 Е larae heat transfer gradient near the upwash side of 
the vortex. The rise in heat transfer rates, although 
demonstrating localized maximums, persists along the length 
of the test surface at a Z value of approximately 2.594 cm. 


When the vortex generator was moved to Ὡξ-δ. 5 


25-32120 Sie 
-6.06 сп, Figures 51 иго ο mE οτε that Ше oS 
overall qualitative data trends are present. However, 
drastic local quantitative changes occur because the vortex 


location is changed with respect to the film cooling holes. 
Local increases in Stanton number ratios are again present 
along with decreases on the upwash side of the vortex. These 
decreased St/Stg values persist not only in the direction of 


flew But also in the spanwise direction. 


Cn 
FJ 


БЕ: о ratios "gre seen im Figures 54 and ΒΒ. 
ο Пе +2 side of the test surface 54/51; ratios are as low 
ав 0.85, The vortex appears to push the coolant from its 
upwash side causing the coolant to disperse in a fairly 
uniform manner. In contrast, without film cooling, decreases 
in local heat transfer near the vortex upwash side are тоге 
localized. 

Figure 59 shows the temperature distribution in the Y, 2 


ӘРЗПЕ БІ Yoertex 3452 with film cooling, Ап approximate 


location of the thermal boundary layer is identified. ΤΠΕ 


results used in this figure area only tentative. 


IV. SUMMBERY ВЕ СОНЕТ 


Baseline measurements show excellent  aareement with 
Stanton number correlations fer a flat plate with constant 
wali heat flux and unheated starting length. Results with 


film cooling show expected trends, and results with an 


ІЛ 


embedded vortex show excellent agreement with data from the 


literature. The effects of the vortex on heat transfer in 
tne Tilm cooled boundary layer are significant and 
impor Sonus ШІ? on the downwash side of the vortex, heat 
tranzefer is augmented, effects of film zooling are negaāted 
and local "Ret-spate" will exist in engines: (2) near the 
upwash side of the vortex, coolant is pushed to the side of 


the vortex, appearing to augment the protection provided by 
film cooling: and (3) as the Vertex location Is change SUENE 
respect 5 film com noles. signifizant ] оса] 
quantitative changes in heat transfer occur even though 
overall qualitative trends remain unchanged. 

Iu ie recommended that in subsequent experimentation 
that spatial resolution be increased, provided it can Бе 
done without compromising the integrity of the heat transfer 
surface. In order to more clearly visualize the interaction 
af the vortex with ulmo ве it would also Бе 


desireable to employ flow visualization. 
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Fhotograph of Test Surface 
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Figure 11, Test &ection Thermocouple Placement 
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Figure 16. Energy Balance Thermocouple Placement 
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Photographs of Wind Tunnel 
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APPENDIX B 


TABLES 


TABLE 1. INJECTION SYSTEM DATA (10041976) 


FLOW RATE MASS FLUX Εἴιι τπ [ος δρα REYNOLDS 
(m^3/s) (Ка/т^ 2+5) (deg C) COEFFICIENT NUMEER 
0.701E-2 0.915Е+1 18.8 0.730 4504.5 
0.654Е-? 0.854 ΕΙ τυ 0.718 4204.2 
0.0ЛЕ-2 0,793Е+ 1 18.8 0.709 3930.9 
О ЕЯ Е ОС 18.8 0.698 оо 
0.51i4E-2 ο TES 18. 8 0,591 3303.3 
Ее UST OE* I 18.9 0.679 3003.0 
О. а21Е-2 0.548Е+1 19.2 0.694 2702.1 
0.374Е-7 0.487Е+1 ` 19.2 0,74 2402.4 
0.327E-2 0.426E-*1 [ove 0.658 м 


TABLE 2. INJECTION SYSTEM DATA (Toe®32 °C) 


ЕО RATE πο LA ΠΕΙ DISCHARGE REYNOLDS 
(m^3/s) (кш ш лс) (deq G) COEFICIENT NUMBER 


ОТЕ 2 


"mJ 
! 


со 


4004. 


O -3 
зо 


0.898Е-1 T Ó 7 
НЕ -2 ВЕ ео 0 9 ЗЕ 
Е OQ.T7T5ET1 το. moat ЗЕ 
ЕЕ 2 Е! 32.0 [D LS DUUM 
ου. -- ΙΙ πη осет E 
Ε/Ε-2 ЕЕ | tro m Ое Е 
NED Quum cim zc i SI DES 2402.4 
ο ευ ο. 4 ΤΕΙ I OS ΙΕ 
Ее ТЕ-2 Ее ος ο οἱ” 1868.5 
С. ЗТАЕ-? 0.477Е-1 στα OS SS 5 
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ο τες ВЕ! ESE D 2659.29 
Cee ie. ge шт! EET COO 2925,3 
а Gn ioe i Sum 0.56988 RE 
EM. OE | нае ΞΕ ЗАТО. 1 
See οι σαν: a 082313 po cu 
ΙΙ Ξ-Ξ ου... CNN DON 4004.9 
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TABLE 3. ΙνΝνΕςστιο μι ΕΙ ' πο ЗЭ ХӘ) 


PEGW erate MASS ee te ENTHE E: DISCHARGE REYNOLDS 
(πι 9/5) (κο πι 2 te) (БЕСТЕ СОЕБЕ ТЕ ВЕТ NUMBER 


ОИЕ Ὁ ΒΡΕ] 43.6 0.538 3603.6 
ο σπα z ο δα 1 43.2 ο 5 ЗЕЕ 
Gg I ELT JG ITE да 2.2651 ое | 
От. О.в νετ. 43.5 OE Pat n. 
O GIRE TE ο νο и 835 0.649 252 0 
О.Е ο 5 4315 0,638 2402.4 
О 2202 Ос ИЕТ 43.2 0.628 Ug eam 
ОЕ О. ODE. 43.4 О.Е λα 5 
FES AS ο πο ЕЕ ques ОЕА" Lol 
U EOE ES 0.450841 εν OQ us s lS 
о. А2 Ег ОЕ Е 43.9 0.433 lom E 
Q που > ο. ο BETI 45.5 0.545 ое 
Оле еа ОПЕЕ ЕЕ ПІ ЕКЕ Doe кз т 
ο ο AS V OTOC E] 43.5 О. 555 2 
ا‎ G:a TERI 43.6 Оа 3d xs 
о ^, S SES 42.5 Qon Eun CU 
ο Е 2 ο σος 43.3 ο. 360355 
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TABLE 4. INJECTION SYSTEM ОАТА (Тс 55 С) 


E COW RATE MASS ЕЕ. EIT ТЕНГЕ DISCHARGE RE οἱ ο” 
(π 9/5) (К9/м^ 25) (DEG ©) COEFFICIENT NUMBER 
ТЕ: — 2 ο πι ΤΕ 1 54.4 oraa: 3276.0 
DAE - ο οσον 343.6 0.542 Есетов 
Oreo ve 2 On (USE 24.1 Ce ee E 
ο E= 2 ο ου πι 54.1 ο ο 8 £620. 5 
ΕΓΞ ЕЕ! 24.4 Oe oss 2402.4 
στον O ο. Ες ο. Dos 2184.0 
Oa cies 0.486Еғі1 EIC 0.606 19597575 
p 374E-2 а ЗЕГЕ 99.4 91 2220 ITA. 
ШӘНІТЕ-2 un toe td απ 0. 581 1528 55 


TABLE J. 


MATERIAL 


INSULATION 
INSULATION 
ЕШЕЗБЕРЕШЕ 
СБЕЙЕЕ 
SIERE 


MATERIAL FROFER TDRS 


ΗΠ ΠΡ, 


0.04 
0.04 
0.184 
OE) 


15.0 


FOR CONDUCTION LOSSES 


$c 
о. 
ο. 
в: 


= 


ea 


п) 


04897 
0843 
416 
Onze 
SES 


соо |Б 
с 


— 
(52 


e 


(т) 


. 0234 
. 0254 
. 1524 
„Ое 
. 00254 


TABLE SS CONDUCTION Εις 


POWER (W) Qup(W) QcondoW) AIC) 
πο 2.07 9.98 3.8 
ШІ” 28 Ber 11.13 13.1 
Eo 3.78 ета с 
ao. a4 ery тает η 
20.45 5.62 14.83 E 


κς Γι ριΡΡΕ ΕΙ) ΑΙ UUVIULHNMEN I EXPENSE 


ВРЕТ АС 
UNGERTAINTY ЭМВЕЗЕТЬ 


Uncertainty analysis was performed using the method 


proposed іп 1953 by KRlime апа МЕСІ таев IRE. 200 ЕЕЕ 


the root sum square method, where the uncertainty, We , of 
some functions is a fonction of the independent variables 
v 
“т . 

1; 


Е у] | C1) 


To calculate the шпевргоуагио MET tne heat transfer 
сое s | Sr ne the to low mu independent variable 


uncertainties where determined: Va Ее W/m p Wry =-0.5 Ge 


ο : : : 
and Wio 20.1 С. Here the uncertainity of convection was 
based ona 2% error in radiation losses and a 1% error in 


conduction losses. The uncertainty of Ty is higher than То 
due to higher uncertainty inherent in our calculation οί 
contact resistance. From these parameters the uncertainty of 
h was determined to be *0.73 U/m^K or approximately 2.5% 
based on an h value of 30 U/m*K. 

To calculate the uncertainty of the Stanton number, the 
uncertainty of the heat MAT. coefficient along with ШЕ 
оТ Таспа independent variable uncertainties were 
determined: W0 9:995 kg/m, Uy. -*0.1 m/s and “ср: 2 ον. 


Ihe uncertainty of Ср was based эп thie assump OT o 
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mens Lan. 


of the 


Properties.: 


From these parameters the uncertainty 


Stanton numbers was calculated to Бе *1.9E-4 or 


approximately οὐ based 


on a Stanton number of 3.6E-3. 


να то в 


= .. 


m МЕ sa ag س‎ z - 
moe ew ewe 8298 ce ace 


.-- io fief OF eS ва” 


ο) 


m 


ا 
be‏ 


ا 
FJ‏ 


рз 


Ν7 


Q1 


а? 


ΕΙ Р 


SOFTWARE 
NOMENCLATURE 
P amb 
Т атр 
Foc Го 


thermocouple temperatures iy) 


Tw - То 


zontact resistance temperature correction 
Vc 

Кед 

hg hf 

spanwise averaged ho: hg 
zurrent 

run number 

q convection (W/m?) 

q comitem Uu» 

q radiation (W) 

power (W) 

ρω 


Юс 


Әсі 


ΤΠ 


--- 


ico) 


S) 


Sit) 


T1 


T7 
v 


ας) 


[i 


Ке, 

Stanton numbers 

tassel ine Stanton numbers 
spanwise averaged Stanton 


То 


voltage 
downstream distance (m) 


contact resistance (K/U) 


numbers 


κ ee κ Ρ ә Ὃν Ὁ ΠΗ ΤΠ 


ж. 


a-s io Фа £ BDH DD теє 


REPRODUCED AT GOVERNMENT EXPENSE 


= 
παν, 


rT 
he 


lO OS) OL fe GJ hac 
m (E am amo ei 


— 
[κι Κι 


f ai pO DO rà w — — — 
СЛ +. (J D — CG D £O —J 
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ί‏ لا ال 
CO Tg‏ 


REM PROGRAM ZTDHTI 

REM 

REM THIS PROGRAM 

REM ACGUIRES MULTIPLE 
REM CHANHEL THERMOCOUPLE 
КЕМ πα 

REM IHO ЕЕЕ Е ИЕ 
РЕМ τοπ ο. 

КЕМ 

БЕМ ІГЗЕНМІУЕОКТІСУ/ЛІСЕРН VER 
SIGH 

FEN NOVEMBER 12z6 

REM 

КЕМ 


ες 119) 


PEM CH UO STS TAAT 
РЕМ COFPER ΕΠΗΞΤΗΗΤΗΗ THERMO 
ГІМІРІ Е> 
REM 
РЕМ 
REM 
FOR Ішіп ТО 139 
} 


КЕМ 

REM EHTER ΜΠΕΣ ΠΡ. DONIS 
КЕМ 

КЕМ 

UISF "ЕНТЕР RUH NÜ <¿HMHMDDYY .H 
ΗΠ 

INPUT Ho 

ПІЗЕ “ЕНТЕК CUREENT CHAMPS 
IHFUT 11. 

ПІР “ENTER VOLTAGECHOLTS)" 
IHFUT V 

DISF "FUH #=",H9 

DISF "ICRAMFS2=", I1 

ΠΙ5 1 ο...» 

КЕМ 

КЕМ 

REM ΞΗΤΕΞ ΠΕΕ Р, AHE Сато 
РЕМ 

КЕМ 

DISP “ENTER FAMECIN HG2" 
INPUT A 

DISF "EHTER ВЕБ ο τι 
РОТ ве 

LISE ENTER AMEE 

IHPUT ni 


124 


осо 


ЛП 
S e c Cu Cad 


f] Gi nie cos 


gov σι ση с 


( 


Cu CU C Ὧι 


Á со іп ca ~j o 


C i σι σ 


2474 74-4 —4 БОРИ? 7 σι 


СП +» сы) [М 


Go im c 


M 
a. 


J 7 74 


Cn Gy hee Caec Jom 
CD Gorn e CA (S Cn emu 


CÓ CO cà CO со с) 


GJ P - @ фи м 
IS Gr iS C i OI GC O 


ош ا 0 ا ترا‎ LO U con 


Со Уа. сл L. 
Ca ο (πι wD 


А — — А — — — — 
μα + aut be pee pee کا‎ 
n TL ¿JO — in O 
Сы D Gc CO Im 


ОТР σα IHN dis = "eA 
UIST ЕН Bona n. P 
DISP "THERMOMETER FRMB<C 3" ят 


ОТЫР "ҮНЇТ τεπητοχί 
РКІМТ "STANT OH Hü RESULTS" 
PRINT "RUN #" HS 


PRIHT 

Lji=n 

REM COHTIHL i 

о. IHUE THE LüüOP 
REM 
OUTP IT НЕ ЕТ. тт] 
ШШЕ ES 

БШЕШЕ 21205 
ет 

ΕΕΕΕ 15.18 
REM 

LEAP 

REM 

FEM 

РЕМ 

REM COHTINUE TO Loi 

Т! ο рр!" (ПОР FOR ALL 
REM 

IF IL S ПЕТЯ 

ГЕ (1-73 ТНЕН = ея 

ВЕ [125251593 ΙΠΕΗ БОТО ;2H8 
GOSUs 2226 

EEEF 156, 1 

Вей =11 

РКІНГТ Sani TC Cos", ACO 
РЕТНТ 

REM TRANSFORM TO SI UHITS 
A=HA¢ S285 Өг 
F =C3233_. fF 
Е ега 15 
Рі=я г 2агжт1? 

= eRe) oo 

РЕМ 

РЕМ 

КЕМ 

REM PRIHT GUT CATA 

РЕМ 


FET JI ΡΕ CKO; 
РЕТНТ Е! 


PRINT "МЕБОСТТУ CM/S)" 
PRINT U 

PRINT "FAME CHIMA?" 
РЕТНТ В 

PRINT "FS TEMP CKO" 
PRIMT TI 


DEI “ТІПЕРМОМЕГЕКЗИШЕСС)" 
DERI πι 


|25 


...%. 


REPRODUCED AT GOVERNMENT EXPENSE 


pM к e БА А н р غور نر نر‎ 2 paad 
F2 FO P2 P2 P2 f PO P POL PO 0 m m n 
D go I VOT fe Го + (Su LO ο ل“‎ 


Fi E MX Sy ا‎ 


(y OI + oJ fà ات مم‎ 


А н — — — — — 
CJ 64 6 Cod Cod Cal Cod 
DNE кк 


Б — — — — 
+ ης 
m T LD Qa =] 


бр 


1420 


p 

“а 

CO-OP hoe CO Uo cO πι 
T Cus 


Cim ccm y Jy QY T C ¿f Cn Cn Cn ¿n cn Cn Cn ¿n ¿n p F f 
Ü J < Q à ii σις δ Ç; (y i Ca O iS: (i (Si ih (Si mnm 


b А н qu je e А KF „ш н — > А -- „л ы БА 02 А к 
Оо ца OT fe lt ои 


тем Tı Ji 


T Er E 
Gi un co 
1%) т ГП 


GI CE Cac < (S τι ΚΟ (5 C cm Cum 


РЕМ 
REM AG PLATE TEMP 


REM 


FOR 1-1 Τῇ (0 
IST C00 pp 


Т "НУС PLATE ТЕШЕ, МЕНЕП 
Кер. суу 
ΕΠΗ TS 

КЕШ 

КЕМ 

FEM ENERGY Επι ας 

КЕМ 


П1-Тз-Ш 

Wl= 9544+) 45401-.89514¢0142+. 
λες 

STAFFED 

Шәшіз2Уу2.15 | 

ОАЕ илипиИиниствьаттр a> π᾿ 
A= | «ῃ 

Q2‏ - 1-1 لا 


БЕГЕН: “ТІ БІНТЕ-І τα Ү" 
РЕТНТ 01 
FIN! “POWER IN CHT" 

INT Ge | 
Τη" “Соне COSS στ... 
ST Gi 
"RAG LOSS <hT)>" 


шг 


HT Ge 
“т "СОНМ 1055 ντ. 
HT W 


“CURR ERTS GR SS 
11 
ο ΤΙ ΥΕ ΗΓ» 


=. 


41414-14 


[ΠΠ 71 21 ὙΠ 2 πι χι а 


жа “Ὁ ΤΙ ΤΙ Τι ΤΙ Ὢ ΤΙ τ ΤΙ ΙΙ TI TI 
— — aX — — -- 


[ м 
КЕМ 
REM 
а= {иил 


Е WALL TEMP. CORRECTIONS 
ETI 


Z-.016 


Sanz 
W=R4. 4897 
КЕМ 
FRIHT "AVG PLATE TEMP, CORRE 
ЕГЕ σου. 
FRINT T323-(15 
РЕІНТ 
ЕКІМТ "Ho 
H" 


(Жи) 
е 
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— I IN On oe C Ira 


DODD D S UR AO O (πι G) (Su 


σι 


00 O CÛ CO ÇO CÛ O ^J ^J ~J ^ ^ ^ JJ 


74 03 CD .E- GJ F9 


i b de^ А au de — — м — — de — — — — 
Cu iS 


— 
сс! 
ea 

MN 


EN 
M 


Iti LO hg CO 
Га на بک‎ Кү 
2 πι σι! 


E E ыл bat 
V 
= 


< 
„М 


[Ó + — н — — — — 
Pm O tO t E E G oa 
MO oo ~1 ON ا‎ С) 


2918 


2128 


(Л ОГО н MW OI oon Js CJ 
بک اک کا ات اس‎ Yi ih Gi i iÇ, CO S: iS; 


Г. 
I Ni NI Ро Мн нн mar 


Го Г Г) Г ald rahy poping ty 


КЕМ CALCULATE НЕНТ ΤΕΗΗΞΕΕΡΕΡ 
СОЕ = FICIEHTS AHO sTAHTON NU 
HEERS 
FOR 151 ΤΗ 2 
IF ісі THEN GOTO 1750 
ΠΠΤΩ 1558 
М-і 
М1=23 
М== 3 
СОТО ага 
М=1 ай 
Мі= 1889 
М=2=2й 
РО I=M TO Mi 
Ci¢cls=ciciri>-O8 
Е Пи e SC {143 
9190140832: 
JisT-H2 
PEINT USING 135980 ; 11,5 (1). 
Hi I 1 
ETE ση ο ED DODE: 2A SDD 
ОПЕ 
НЕХТ i 
NEST 4 
РЕМ 
КЕП ТРЕЕ ШЕЕ ЕЕЕ СЕТМПІ ГІ 
2 HUMSER 
REM 
FRIT 
КЕМ 
Hl= ΩΗΗΕΊ1ΞΙΕ 
Εξ] “Η1 
Kic 1151 113289 
=. СВЕН 
2163 1=1. Заз ая 
21(32=1. 459669 
м1023=1. 6999 
51652-1.7651439 
AICS 91355 
“~1¢Sas2¢ .466623 
РЕІЧТ "EOW 8 РЕ" 
ΕΡΕ Ісі Tu 3 
ЕРСІзсХіІ СІОЖҒ 
PRINT USING 2166 ; Ι,ΒΖ(Ι} 
IMAGE DD, 3⁄4, D. 4DE 


НЕХГ I 

| 

I1 CALCLILATE AYG 
GOÜSIIE 2426 

ENC 


i 
SUBROUTINE 


STANTON NG 


PGE TO TEHPEPRITIJ EE 
Е. 


ο ο πα πο 


ео 


=< + PD 


12] 


а 44! 


= 


GOVERNMENT EXPENSE 


| 


n3 
[x 
σι 
δι 


[.1[. [2 
fe ro ro 
Uu Οδ πι) 
ا‎ 


S280 
221686 


ا ت 
ΞΡ‏ 
2241 
2358 


b q 
ει σα اٿ‎ 


T 
=$ 


+. pa Kal Ci taj Cd 
1. ym بآ‎ 
I 


Pato hs Pah ti 


г 
— w 
m Z! 


2428 


τος РЕЖЕ 1.ЭЗЕ ЕБ ӘУ ЖЕЖЕ+. ЯЯ 
CSI EEES 2l TE ETETE 
RETURN 

l 


| FIHO THE AMBIENT TEMFPERAT 
URE 
| 


пит μι οι πι ο ΤΙ 
и 


ΕΤΕΠ οσο 
ine ΤΕ 217A 
U=T 


GUTPUT Cle > “Hl sti aa Vit 


т Ul 5 AMBIENT TENT. 

! 
ЕЕТМЕН 
КЕМ [ΠΞΗΠΕΤΕΗΙ TENE, 
КЕМ 
КЕТЫРН 

| 


! SUBRRDUTIHE TO FIMD Βυῦ ST 
ANT OH На 


FOR 1=1 Tü 7 

НІС Із-й 

M-I£1i-198 

FOR !=М ТЇП [Ж11 

HII szEDCIoTESJS 

NEKT J 

О το 11 

Бїс ΕΠ σπα το 
НЕХТ í 
Н1 С 2 2 =8@ 

FUR i-?2 TO 7 
H1z35»-H148»5-HÁID 
HEXT 
FOR 1=1089 TO 188 
HiC(2»5-H1(CB8»5-4HXXID 

НЕХТ I 

Η1 3 ΞΗ1Ι (Βλ 11 

оза ен P αν 

| F2SIHT RVG ΞΤΗΗΤΟΗ ΝΟ 
PRIHT 

РКІМТ ge 5 AVG ST ΝΟ" 
E 

PRIAT Us ING "2696 зо 
IMAGE 00,35, 0. 40Е 

НЕФТ 1 

КЕТЈЕН 


es WD 


id 
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428 
438 


535 


556 


REM 5ΡΠΩΞΠΒΗ STOAT 

REM 

REM THIS PROGGEHHM 

REM ACQUIRES MULTIPLE 

REM CHHHHEL THERMOCOLIPLE 
REM DATA 

REM AHO CALCULATES HEAT 
REM TRANSFER 

REM 

ο οκ ы ЕЕН МЕР 
Так 

РЕМ ЧП ЕМВЕЮР 13&6 

КЕМ 

КЕМ 

КЕМ 

CLEAR 

DIM 1.1349). 158) 

DIM Βί 1.» Ες 115) K 118) 
DIM He ος a 581187) 


£7 598 


REM ENTER AMPS AND VOLTS 


REM 

DISP “ENTER RUN NO CNHDOYY.H 
HIM» " 

INPUT НЯ 


OISF “ARE YOU HORKING (1 2HIT 
Н ПЕ (29 WITHOUT FILMCOOLING 


IHPUT 21 

IF Тісі THEN GOSUB 2354 
DISP "ENTER CURRENTCAMPS >" 
INPUT I1 

015Р_"ЕНТЕК VOLTAGECVOLTS>" 
INPUT V 

DISF "KUN Hc M3 

DISP "IXRMPS2-z", I1 

DISP "/CVOLTS5z",V 

КЕМ 

REM 

REM EHTER GEL P- ANB CNDTNS 
РЕМ 

КЕМ 
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GJ PO — Cpu Ig,‏ > اک 


J N D Tı J, Ti Ci Os Jy CT CR Cn СЛ 


сос GO 0 O ο c C m © 3 © © © 


b Og os‏ مہ زم اہ ہچ 


e 


=I “INI М = 
Cn n 
бо 


ISF “ENTER PANECIH НЫ)” 
HFUI T 
ISF E d DEL РА Haw 
INFUT 

ISF "ENTER ТАЛЕ СЕК 
INFUT яі 

DISF " PRMB« ІН Нб>-".Н 
DISP "DEL FOI EUS E 

DISF "THERMOMETER ΠΕΡ, Π! 
КЕМ 

КЕМ 

КЕМ 

КЕМ 

ВЕ "НЕЕ СОН) А 

PAUSE 

PRIHT “STANTOH NO RESULTS" 
РЕІНТ 

PRINT "КОН #".НЭ 

PRIHT 

1=й > 

КЕМ СОНТТНИЕ ΤΗΕ LOOP 

ШЕРТ 


OUTPUT 7/09 ;j"RI“%“;jË De Pl ΤΙ 


C AD 
Τι 


D 
( 


0 


BEEP 18.18 

КЕМ 

CLEAR 

REM 

REM 

КЕМ 

REM CONTINUE TO LOOP FOR ALL 
THE2NOCOUPLES 


REM | 
IF 11<79 THEN GOTO 760 
IF L1=79 THEN 15 τα 
ІҒ 1152139 THEH GOTO РЕВ 
ООБИВ 2416 

BEEF 158,198 

НӨ»? =! 


PRINT "АМВТЕНТ ТСС>:", Вс 
PRINT р Ge? 
REN TRANSFORM TO SI UNITS 


ΗΞΗΣ 2455 


ь 


ΓΞ Κ548.Υ 
Тіс-Сі(16924%223.15 
Еі-В/с28ҮЖТ12 
U1zi2tC^R125^.5 

IF 71-1 THEM GO0SUB 2940 
REM 

REM 

REM 


REN PRINT OUT CATA 
PRINT "“QENSITY εκ σσ" 


PRINT R1 
PRINT “VELOCITY ch o 
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=J OVN f 
i Cn G (2 


б> mA m ны на А pt А А 
à P2 F3 P2 PO P) f а в ee مر‎ 
(e X n ( (СД < 


ы Ізі 


a» ub TJ UC 4 G1 P9 — m LD οὗ 
АА 


CJ 6J Г) РГО ГО 
tel lo ©0 ا‎ 


هلر مغل م خم عارع مرم mM‏ 2 هر غ خم ملم 
. 


P aJl ta ل ف‎ aje 


Со CO ay (Ji Cn T+ 
TORCA m m s (& C) C C um 


к. 
+ 
م‎ 


1428 
145 
1440 
1458 
14615 
1474 
1430 
1+ 
1544 
1515 
1528 
1538 
1548 
1558 
1568 
1578 
1588 
1598 
1648 
1618 
1624 
1638 
1644 
1658 
1668 
1678 
1680 
1690 
1700 
17168 
1229 


PRINT Ul 
FRIMT “PAME έΝ/ Μο)" 
PRIAT 
PRINT "FS TEMP «K»" 
PRIHT Ti 


FRINT "THERMOMETER AME€C >” 
PRIMT H1 

REM 

REM AWG PLATE TEMP 

REM 

REM 


= τι 


PRI, 47 "AVG ΕΙ ΗΤΕ TEHPXC»" 
FRIAT T3 
КЕМ 

КЕМ 

REM ENERGY 


REM 


EALANCE 


1.45*091-.@51+%01^2+. 


clo9KCRS^d-US9^4) 


FRInMT 
FRIHT 
PRINT 
PEINT 
РЕІМЧТ 
PRINT 
PRINT 
РРІЧЯТ 
РЕ I ЧТ 
FREIHT Q 

FRIGAT "CURRENT CAMPS)" 
PRINT 11 

PRIHT "VOLTAGE 
PRINT VM 
КЕМ 

РЕМ 

КЕМ 
M-LSELIUN. 

REM 

REM WALL TEMP. 
РЕМ 

2= .315 

05=9%2 

(=G 42897 


ΓΗ ΤΗΤτε- HPB C5" 
01 


"POWER IN CWT" 
03 

"COND LOSS ЫТ" 
ῃι 1 

"RAD LOSS CWT)" 
Q 2 

"CONV LOSS CNT)" 


CVOLTS>" 


CORRECTIONS 
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кекке M = خسم خسم‎ 
СО о со оо со со У ل“‎ 
+ Ы г = (c b 03 

J O GE S CO (IS iS επι 


j 
с 
a 
х 


сс 
сл 
c 


Ξ86Η 
2074 
2034 
20894 
216% 
2118 
2128 
2130 
2148 
2158 
2166 
2178 
21880 
2198 
2200 
2218 
2228 
2238 
22408 
2250 
2688 


КЕМ 


ЕЕІМТ 
FRIST “Wo Η 

T/STH" 

REM CALCULATE HEAT TRANSFER 
COEFFICIENTS ANT STANTON HU 
HEEFS 

| OPEN FILE TO RETRIEVE BHS 
ELIME ST VRHLUES 

RSSIGN# 1 TQ “HOATA" 

FOR J=1 TO 3 

IF 42-1 THEN GOTO 1826 

ΠΌΤΩ IRER 

i= j 

MN1= 73 

Ме= и 

ООТО 1590 

М=1ійӣ 


i 


ж 
к= 

| 

حا ل 
UR.‏ 
ол‏ 


πς dG O DEF CT 
zg ^ Mr O RN 
i مم ر“ سم ,حر‎ «SW || | 
Amr: 

c 

os. 

(=a 

ы? 


liZIHG 1928 ; J1.,H«I». 


AJ 
4 
ut 
ed 


seco 


IMAGE 2D,2xX.S2D.UDDE, 23, DID. 


ар INT# 2 2, у =. 
SIGH# 1 TO % 
ASSIGNS 2 TÜ Ж 
КЕМ 
REM CALCULATE LOCAL FEEYHDLD 


S HUHBEE 


КЕМ 

ΕΚΤ 

REM 

Η1ξ. ΜΗΒΗ 196 
F=U l 1 
Хісісзгі, 1559 
«289563 
. 39429 
. 43669 
41t5?-1. 60989 
А1663-1.76149 
А1<721.913989 
416222. 06629 
PRINT “ROW # 
ҒОК 1-1 TO 8 
ЕгСІозтсхі СІОЖҒ 
PRINT USING 2230; 
IMAGE 00,3%,0.4DE 
НЕХТ 1 

| 

! CALCULATE AVG STRHTON NO 


КЕ” 


1,Ё2‹1› 
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Γι) ΓΥΡΟ 
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CJ CJ 56 CI po PL PR 
Ü] [O — tS LO (9 ~) 
(τι (τ (πι (Sh Q iS ὅτι 


2340 


λος) 


(n ¿n Cn n (n р. 


gicn 4» oro cu 
C < i {Ὁ اک ات‎ 


Го го Го Го Гого Гага 


T سک اک نت لک‎ C 


Ρο ΡΟ ΡΟ ΤΟ ΤΟ Fo 
(Тэ Сг КТ emo (n Cn on 
fe ph yo co mJ 


2098 
283969 
59518 


FRINT 


сае ЧЕЗ 


ΕΗΏ 
i 
SUBRÜLITIHE 


| 
і 
l 
i VOLTAGE TO TEMPERATURE 
|! CONVERSION 
| 
| WEMY Tü Ter) 
l 
Εξ} πα 
Т-25 5ҒЗЖЕ-| ΞΠΡΟΡΑΣΕΧΕ. 99 
= ЗУБ ЧЗ ETETETE 
RETURN 
| 
| FING THE AMBIENT TEMPERAT 
RE 
| 
OUTPLIT TAS ;"RI";Ci£B»;"UTIl 
ЕНТЕР таа s 
GüzlBg 22йй 
U=T | 
UUTEUT 7O3 :"AI";T1CBD: UTI] 
I il = AMBIENT ТЕМЕ. 
RETURN 
РЕЙ X II-HHBIEHT TEMP. 
REM 
RETURH 
i 
| SUBROUTINE TO FIHO AVG ST 
антон NO 
| 
FOR 121 TO F 
Hl Із =й 
М= 1511-10 
FOR М=М ТП 1Ж11 
HICISSHiCIO*HCJ) 
HEXT J 
HICIO?-HiCIS711 
SICI-HICIDOAZCR1XULXCO9) 
НЕХТ 1 
Η1 823 
FOR I=78 To 73 
H1C332H1 CS» *HCID 
HEXT 1 
FOR I=146 TO 148 
H1C3»-2H1C8»*HCIO 
HEX 
H1C82-H1C82711 
21C332HL1CSO/CRIXU1XC39) 
| 
| PRIHT AVG STANTON NO 
РЕТНТ 
PRINT “ROW # AVG ST HO" 
FOR IT=1 TQ 2 
UST 20 TCI 
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O OI ON On‏ لت حسم 
NO ei co S‏ 


Po Pl m mA mama 


2208 
3260 
ЗЕ 

3280 
3290 
3300 
5310 
3326 
5530 
3360 


INAGE DD.34,D.4DE 


HENT I 
КЕТПЕН 
i 


| SUBROUT DHE FOR FC ΠΡ 
DISS "ENTER FLU VALUE CAD 
[МРАТ Fi 

DISP “ENTER FEEN UEC έτη 
Hzn: 

INFUT Fi 

Fi=Pit248.7 

Е} ==1%. 989992455 

ЕЕТМЕН 

| 

[I SuERUUTIHE FOR ο CALCULT 
ЕЕ 

ТС ΠΠ» 
23 

! LOHVERT PLEH TEMP TO IWJE 
ΕΤ ТЕПЕ 

Te-i 45463xT5^.9695162 


Τοξτε ο 3 i 
4 


RS=A’CTY 287?) 
Us3ZiFIXZÁ/^R32^.5 

03ЕРЗХИЗ 

04=02 072 | 
F2-u2€t.0609525^» 8006156 

i PRINT RESULTS 

ΡΚΙΗΤ “INJECTION TEMPCC>" 
ΕΚΙΗΤ ТБ 

PRIHT “PLENUM TEMP«C»" 
PRINT TS 

FRINT "PLENUM DEL P «М/Мео> 


F1 
PRIHT "KHOC «4KG^M^25" 
R2 
VT "VELOCITY,UC (M^$)* 
! Е, 
NT "MASS FLUX (KG/M~24S >" 
NT “REYNOLDS NO. DIR)" 
E 
PRINT "DISCHARGE CUEFFICIEN 


PRIT O4 

PRIHT "BLOWING RATIO" 

PRIMT D2^7cR1XU125 

FRINT “DENSITY RATIO" 

PRIMT KZ^/R1 

PRINT "VELOCITY RATIO" 
PRINT U2^7U1 

PRINT "MOMENTUM FLUX PRAT IN" 
PRINT K2XU2XU2^;COR1XULTIXUI125 


RETURH 
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259 


268 
гей 
234 
290 
380 

218 
320 
338 
346 
258 
η. 
aff 
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REM PRUGRAM TPROF 

REM 

REM THIS PRUGRAM 

REM ALGUIRES TEMPERATURE 
REM PROFILE OATA USIHG 
PEM A THERMOCOUPLE 

REM МЕ 2418 

REM 

REM LIGKEHNI, HOL 1556 
REM 

CLEAR 

О1М Сел, ТС: 
DIM ASAD, Ες 5Η 
WM ot 5h yet ae 
НЕЙ! e M 

REM 

REM 

[CIE IS 

REM 

REM CHHHHEL 1365 FOR TC 
КЕМ 

REM CU-CONSTAHTAN TC 


Уи 


FEN 
REM 
WISF 
ΗΠ" 
THPUT No 
HISP “a EUERT LON si) > © 
INFUT 41 
ПІСР "2 LOCATION Мә" 
IHFUT 21 
DISF "UIHF M/SEC)" 
INFUT Ui 
τ tee лды 
ΠΠ. C2 
ΠΙ5Ε "à CONV" 
INPUT Wl 
DISF "THRLL <Ü)" 
INFUT ТІ 
DISF "OEHSITY εκ Πο)" 
IHFUT 01 
П15Р “СР (WHT SECZKG K>" 
INPUT Ci 
GISP “WALL DSPLMNT <ММ2" 
INPUT Y3 
 . 

=Ü! EIC. 23 
т Sao” 0156112) 
μις ΒΜΗΜΡ156 
REM 
REN 


= 


'EHTER RUH ND CMMNODYY.H 


REM 
PRINTI "TEME FROFILE STUDY HO 
ΙΕ. 
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id 
IHT О "ЖЕХЖЖЖЕЖКЖЖЖЖЖЕЕЖЖКЕЖ 


p5-50.3 
REM 

REM 

REM ENTER MAIN LOOP 
КЕЙ 

FOR [=1 ТО 99 

КЕМ 

Sr Ее 

DISP “ENTER ποτητ ο - 
TREUT ЕЕ 

DISP "ENTER ПЕСЕББЕ» 
INPUT RZ 

ΤῈ 
y2z-Yi-mzt25.47i260*$132 
IF I-1 THEN οι 5 

IF I=} THEN БОТО 799 
ете еә 
КЕМ 
OTs 


Е 
F 
E 
ISP "HIT CONT. FOR D.A.“ 


3310 


Е 
ЕБ 
SP “CHANNEL #"5C1¢L1) 


FOR K=1 ТО 29 
& OUTPUT 789 ;"RI";C1C(L15;"VTI 


Wk s 
NEXT K 

X= 20 

PRINT "EU D; σου, 
PEINT USING ^D.6D^ ; K 
GOSUB 2220 

ВЕЕР 10.18 

CLER? 


J REM 


fe je‏ — م نر نر نر ن ضا 
p e ob (aua Pe i (La je (wx‏ 
cn fe OI AD >‏ 03 4^ 03 
Сосо Q @ gs‏ 


а а 
nro n — 
m cu 
(Ὁ εξι Gs 


КЕМ 

КЕМ 

T<I:x=T 

F:E H 

FRINT “POSITION DATA"; 


FREIHT 

FRIHT "EOTRHTIONS";R1; 
РЕІНТ " OEGEEES";R2; 
FRINT 

FRINT "Yz";YCID;" MH"; 
РЕІНТ 

PRINT "Ts"TCIOS ЕЕ"; 
PRINT 

FRIMT "TXRTTTTIELTEPEEDEEERTPEPLTEX 
ЖЖЖЖ" 

FR IT 

КЕН 

ΕΕΕΡ 1, 18 
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1668 
16198 
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1630 
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1650 
1668 
1679 
163H 
1e9n 
1/64 
1718 
17289 
1730 
1745 
1750 
1°64 
1°76 
1280 
1738 


e c CO CS (C 


REN 
REN 
DISP 
BIS 
INPUT 
te = 
КЕМ 
NEXT 
REM 
REM 
r=! 


es 
fair 


CT 


c= 


REM 
REM 
GOSLB 
КЕМ 
КЕМ 
КЕМ 
КЕМ 


FRIAT 
PRINT 
РКІНТ 
РЕТЧТ 
PRINT 
Pop Til T 
РЕІМТ 
РЕІНТ 
PRINT 
PRIHT 
PRINT 
FRIHT 
PRIT 
PRINT 
FEM 
РЕМ 
КЕМ 
REM 
FOR 


“LAST POINT?" 
"1 IF YES" 
J 
1 THEN GOTO 131A 
l 
I 
Io*TC25 o 
2440 
КЕНЕ FROFILE STUDY WN 
5" 
D I а 
ж.ж. 
"RUN NUMBER" 
М 
x LOCATION MY" 
"С OCG <Ha” 
ΠΟ ΤΗΕ «М/5ЕС;" 
U1 
"СЕР" 
p 
"ο Γον (ΜΗΤΙΣ.ΜΖ}" 
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Ες 11η 
Π1 i 
“CP iNT Έτ τοι." 
C1 
"HALL DSPLUCMHT CMM? " 
Y 
"MRBLL TENE DL +" 
T1 
"FS TEMP Co" 
T2 
“FRICTION TEMP CKD" 
TZ 
“THERMAL DEL ММ) * 
05 
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itt IT 
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YCI»^7DS5 
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РЕМ 
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KAK RE RAR RAK 

ӘКІМІ 1 Y CMM» То 
Co” 

FOR I=1 Tū IF 

ГЕНТ USING TS PE СИТ ST 
CI» 


IMHGE ор, 3х, 000. 10Е, 284, 080.2 
(ШЕ 


ЧЕАТ 1 

FRIHT 

КЕМ 

КЕМ 

PRINT 

PRIMI “I Y7 BEE СШНЕ-Т 
2/0] - TRO " 

РОК I=] TI IF 

РЕІЧТ USING 12568 ; I.D«4IO,H 
C13 

IMAGE 00,35, 0. 50Е, 2м, 00. З0Е 
NENT 1 | 


FRINT " I АСЕ T- 
осо" 

ак 1=1 ТО IF 

PRIHT USING 26846 ; I,U«4I»5,Z 
CN 

IMAGE СО, 34, 0. ЗОЕ, 2х, 00. 20Е 
NEST I 

РК МТ 

КЕМ 

РЕМ 

RET 

РЕМ 

КЕМ 

КЕМ 

FRIHMT " [ Y T+ 


ҒОК ігі TO I? 

FRIHT USIHG 2168 ; І,ССІО,Б 
(I) 

IMAGE DD,35,DDO.1DE.2X,DDD. 
10E 


НЕХТ 1 

PRINT 

КЕМ 

КЕМ 

КЕМ 

КЕМ 

КЕМ 

PRIHT 

ОТОР "КИН GO MP ie ies 
PRIHT "EUN COMPLETE" 
РКІНТ 

КЕМ 
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SUBROUTINE 


VOLTAGE TO TEMPERATURE 
CONWERS ION 


l 
l 
i 
| 
і 
| 
| 
l 
| 
l 
F = x 
Tz: 
КЕТПЕН 
КЕМ 
КЕМ 
КЕМ 
КЕМ 
КЕМ 
КЕМ 
REM 


КЕМ 
S 


t 


= kz1 TS 

Ел а -TCKI3/7CT1-T2> 
ς18-ΚΟΞΥΞΕΥ2-ΥύζΚῚ 

ы. K 

FOR Е=1 TO 17 

IF cuK>>Z2Z THEN GOTO 2698 

NEXT К 


REM 
FdsGüiK-20^2XGCK-120*GCOK 2 ^2XG 
(K-25»*GXK-12^2XGCOKD0 -GCK5^2X 
GOK-13-GOK—-19*2*6¢€K-29-GOK- 
=. 23: ú E 5 

РЕМ 

REM 

F1-2'CECK-13-ECK-25D5 KCGCE»-G 
(Κ- 232-CECKD- οκ ОКР 
20-6 i:K-222»-F4 

REM 
SZOCECKO-ECK-2»25£CGCK-1»5^2 
-GQX4-22^2»-CECK- 1»? -ECK- 25 


CC 27 2-GhuK-22^22»5»7F4 

КЕМ 

Faza. k-23^2*6«4K-1»5XESCKO-*-GXUXK 
223%1:С<К-2>2ЖЕСК-192 

F32r3-G(OK-12^2X3GCKO XECK-Z2- 
GCKO^ZXGCOKCC 12 ECK-2) 

Ғасг3-СҰК-19е2% E - 

το ο τ 
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КЕМ 

O5=i,°k—-1) 

FUR k1-1 TO 18 
S-2ü05-Cz-F1xX*D5^2-F2t*DS-F3»-» 
CczXritD5-F2» 
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REM 
RETIRE 
REM 

REM 

REM 

REM 

FEM 
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Бо ae 


F 156: 125901 


11586.1201 


APPEND IT. 


ESSE DNE БОЛ 


NOMENCLATURE 
hens 


thermocouple position 


PPE SOP Ir I ION 


10 m/s 
10 m/s 
10 m/s 
15 m/s 
20 m/s 
temp. profile, OL mMm/ Se A70 03 т, 


= ‚85 т 


temp. profile, Do m/e, Z-0.038 m. 
A-1.85 m 

temp. profile, 10 m/s, Z=Ū m. 
A-71.85 m 

temp. ororile, ШІ m/s, 2-0 п 
Х=1.443 т 
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